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ABSTRACT 



Alms. We report the outcome of a deep multi-wavelength study of the IC 2391 young open cluster We aim at uncovering new low-mass and 
sub-stellar members of the cluster and identifying new debris disk objects. 

Methods. A 30x30 square arcmin area in IC 2391 was observed using the wide-field imager at the ESO 2.2m telescope. The completeness 
limits of the photometry at 3it level are K=24.7, Rc=23.1 and /c=23.0, faint enough to reveal sub-stellar members down to ~0.03 Mq. Our 
membership criteria are based on the use of our optical data, in combination with JHKg magnitudes from the 2MASS catalog. We also 
estimate the physical parameters of the selected candidates. Debris disk candidates are identified on the basis of their infrared excess emission 
using near- and mid-infrared photometry from the Spitzer Space Telescope. 

Results. Our optical survey, which has a limiting magnitude at 3ct level 1-2 mag fainter than previous optical surveys conducted in IC 2391, 
revealed 29 new low-mass member candidates of the cluster We estimate the contamination to be at least ~50%. We constrain the fraction of 
sub-stellar objects in the range 8-15% and discuss possible explanations for the deficit of brown dwarfs in this cluster We also identified 10 
candidates in the cluster showing IR excess emission consistent with the presence of debris disks. 



Key words, stars: low-mass, brown dwarfs - 
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1. Introduction 



stars: formation - stars: pre-main sequence - stars: circumstellar matter - ISM: individual objects: 



Over the last decade the general form of the stellar IMF 
has been observationally established, while the situation re - 
mains less solid across the hydrogen-burning limit ('Reid'2005'). 
More recently, infrared ( IR) survey data ( e.g. Bouvier et al 



the clarification of the dom inant brown dwarf ( BP) formation 

mechanism in star clusters ( Bouvier et alj|200"5 ). 

Because of its proximity ("~1 50 pc: iForbes et al. 2001 ) 



1991 lLodieuetal.l l2005l r2 006a b d Eool iLuhmanI l2004b ; 

Qoet 



relatively young age (30-50 Mvr: iBarrado v Navascues et al 
1999'. 12004 feasri et al.lll99l lMermilliodlll98lt IStauffer et al 



Luhman et al.ll2006l:lMoraux et all 120031: iMuench et all 12003 



2007), have produced rehable statistics for substellar objects. 
These overall results, in combination with optical photomet- 
ric and spectroscopic follow-up observations, allow to quanti- 
tatively address such issues as the continuity of the IMF across 
the stellar/substellar boundary, its invariance or dependency on 
local conditions, etc. Such studies also hold important keys for 



,199&) and i nsignificant extinction along the line of sight (E(B- 
V)^0.006: IPatten & Simonlll996l) . the IC 2391 young open 
cluster is long since considered an ideal location for star for- 
mation studies. Consequently, IC 2391 is one of the best stud- 
ied clusters. There are about 180 members known in this clus- 



ter, with accura t e proper motion d eterminations (IDoddl 12004 



Monet alJl2003t iPlatais et al.l 12007 1. among which a dozen of 



substellar objects. The only determination of the IC 2391 mass 
spectrum in the very low-mass re gime has been presented by 



Spezzi, e-mail: 



Send offprint requests to: L 
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* Based on observations carried out at the ESO/La Silla, Chile under 
proposal number 68.D-0541(A). 



Barrado v Navascues et al.l (120041) . According to these authors, 
the mass spectrum between 0.5 M© and the substellar limit 
(~0.07 M0) follows a power law with a-index ~0.96. Below 
this limit, there is a sudden drop that the authors attribute to a 
local drop in the shape of the luminosity-mass relation at spec- 
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tral types M7-M8 and, partially, to the incompleteness of their 
spectroscopic follow-up below 7^=18.5 mag (i.e. 0.05 Mq for 
cluster members). 

However, other young clusters (e.g. a Per and NGC 2547) 
show a lack of substellar members more or less with 
the same masses as in the case of IC 239 1 (see 
Barrado v Navascues et alJ 120021 : Ijeffries et alJ |2004j). This 



deficit may extend to lower masses or may only be a dip, 
perhaps caused by an imperfect understanding of the mass- 
luminosity relationship at temperatures around 2400-2500 K 
jPobbie et aljf^j02i) . However, the incompleteness of the sur- 
veys conducted in these clusters leaves the question open. In 
particular, the number of substellar objects discovered so far 
in IC 2391 is too low to draw any conclusion on a statistically 
significant basis. 

In this paper, we present the results of a deep opti- 
cal, wide-field imaging survey of IC 2391, complemented by 



IR photometry from the 2MASS catalog (ICutri et all 120031) 
and the Spitzer satellite. This survey, which has a 3cr pho- 
tometric limit 1-2 magnitudes fainter than previous optical 
photometric surveys conducte d in IC 239 1 (i.e. Ic <21- 



22 m ag; Barrado v Navascues et al. 200 Ibt Patten & SimonI 
19961: IPatten & Pavlovskilll999h . aims at uncovering new low- 
mass and BD member candidates, to be confirmed with fu- 
ture follow-up spectroscopy, in order to provide further in- 
sights into the substellar mass spectrum in this cluster. In ad- 
dition, the age and distance of IC 2391 offer a unique com- 
bination to study the evolution of debris disks around low- 
mass stars, which is a crucial step toward our uriderstanding 



of the terrestrial planet formation mechanism (e.g.|Cu rrie et aL 



2001 



20071: iDecin et al.ll2003l;[Dominik & Decinll2003l:lHabing et al 



The study of circumstellar disks has recently under- 
gone a substantial improvement thanks to the data from the 
Spitzer Space Telescope. These data are very adequate for 
disk investigations because the Spitzer sensitivity and wave- 
length coverage allow studies of statistically significant sam- 
ples and probe the in ner planet-forming region of disks (see, 
e.g.. lLada etaDl2006h . Spitzer observations have shown that 
debris disks are found around solar-like stars at a wide range 
of possible distances (1-50 AU) from the central star and tem- 
peratures (10-650K ), with an age-dependent frequency (see 
Siegler et alj|2007l and reference therein). A pioneer study of 
the debris disk population in I C 2391 has been recently con- 
ducted by ISiegler et al.l (l2007h using the Multiband Imaging 
Photometer for Spitzer ( MIPS). Using our o ptical photome- 
try, the MIPS catalog by ISiegler et all ( |2007|) and additional 
IR photometry from 2MASS and the InfraRed Array Camera 
(IRAC) on Spitzer, we search for new debris disk candidates 
in IC 2391. These studies ultimately aim at measuring the fre- 
quency of debris disks around very low-mass stars, in order to 
assess whether they experience a planetesimal ph ase in their 



evolution, as observed for more massive stars (e.g. Chen et al. 


2005: Kessler-Silacci et al. 


I2OO6; Lada et alj 


20061 Riaz et al. 


2006:.Siciha-Aguilar et alJ 


2007: Young etal. 


2004h. 



In Sect. |2] we describe the WFI survey in IC 2391 and the 
data reduction and catalog extraction procedures. In Sect.[3]and 
Sect.|4]we describe our membership criteria, present the sample 
of new cluster member candidates and estimate the contamina- 



tion level of this sample. The fraction of sub-stellar objects in 
IC 2391 is discussed in Sect.|5l while in Sect. |6] we investigate 
the debris disk population in this cluster Our conclusions are 
drawn in Sect.|2l 

2. Observations and data reduction 

2.1. The WFI observations 

The optical photometric observations were performed in ser- 
vice mode using the Wide Field Imaging (WFI) mosaic camera 
attached to the ESO 2.2m telescope at La Silla (Chile), from 
November 2001 to March 2002 (Program ID: 68.D-0541(A), 
PI: Isabella Pagano). WFI is an 8fc x 8fc CCD mosaic, covering 
a 30x30 square arcmin field with a pixel scale of 0.238 arc- 
sec/pix. These observations are also intended to investigate the 
variability of young stars at the bottom of the main sequence 
due to their magnetic activity, which will be presented in a fu- 
ture paper (Pagano et al., in preparation). For this reason, our 
observing strategy consisted of repeated exposures of the same 
field. In order to avoid photometric contamination from the 
bright AB type stars in the core of IC 2391, our WFI pointing 
is centered at R.A.=08'':42"':3P and Dec.=-53'':02"':58', i.e. 
about 0.25 deg east from the cluster center (Fig.[T]), and covers 
about 60% of the sky-area spanned by IC 239 1 . The region was 
observed in the V, R and / WFI broad-band filters for nine non- 
consecutive nights, with 1-7 observations per night per filter 
The seeing conditions were always better than 1 .2 arcsec. This 
strategy ensured a good phase coverage for objects with rota- 
tional periods between 0.2 and 10 days, i.e. the range expected 
for very low-mass stars at the age of IC 2391 (30-50 Myr). In 
order to cover the gaps between the WFI CCDs, each obser- 
vation in each filter was split into three individual exposures 
(ditherings), shifting the telescope pointing by ~1 arcmin be- 
tween consecutive exposures. Since we are interested both in 
identifying the lowest mass cluster members and recovering 
the early type ones, we performed two series of observations 
per night. The total exposure time was 270 s in the V-band and 
75 s in the RI bands for the "long-time" series and 30 s in the 
V-band, 15 s in the /?-band and 21 s in the /-band for the "short- 
time" one. In this way, the photometry for early type mem- 
bers saturated in the long-time exposures is recovered from the 
"short-time" ones. A summary of the observations is reported 
in Table [U 

2.1 .1 . Pre-reduction 

The raw images were processed using the IRAlQ mscred pack- 
age and a number of scripts ad hoc developed both under IRAF 
and under IDL0. We followed the standard steps for mosaic 
CCD data reduction, i.e. overscan, bias and flat-fielding cor- 
rection. For each night in which observations for our program 
were performed, bias and twilight flat frames were combined to 
obtain the night master bias and flat, respectively, which were 



' IRAF is distributed by NOAO, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under 
contract to the National Science Foundation. 

^ Interactive Data Language. 




Fig. 1. R-hand image of IC 2391 from the STScI Digitised Sky Survey. The solid line display the total sky-area surveyed with 
the WFl CCD mosaic camera at the ESO 2.2m telescope while the dashed lines trace the area observed with IRAC and MIPS 
at Spitzer (see Sect.|6]l. The arrow indicates the cluster center The filled squares represent the new member candidates identified 
in this work, while the circled symbols indicate the objects showing IR excess emission (the stars mark the confirmed cluster 
members, while squares and crosses mark the member candidates and the four sources of dubious nature discussed in Sect. l6.1.ll 
respectively). 



then used to correct the science images. Possible residual ef- 
fects due to vignetting and sky concentration caused by the 
geometrical distortion were reduced by subtracting from each 
image a background's 2D-polynomial fit, obtained by using the 
imsurfit routine under IRAF. Fringing strongly affects the WFl 
observations taken at near-lR wavelengths. In order to remove 
it from our /-band images, we subtracted from each /-band sci- 
ence frame the fringing pattern frame available as part of the 
WFl standard calibration plard, scaled by a specific factor to 
account for the amplitude of the fringes in the individual sci- 
ence frames. 



2.1 .2. Astrometry and co-addition of images 



The astrometric calibration of each exposure was performed 
using the msczero a nd msccmatch packages under I RAF; the 



GSC (ISTSc]|l2006i) and TYCHO (IHog et al.lll998l) catalogs 



were used as reference and an absolute astrometric precision of 
0.6 arcsec was achieved, as confirmed by a cross-check with the 
2MASS point source catalog (see Sect. l2.1.4] i. Individual expo- 
sures in each dithering set were then resampled, flux-scaled and 
combined into a single mosaic by using the mscimage, msci- 
match and imcombine tasks under IRAF, respectively. The fi- 
nal stacked image is a 8A: x %k frame where each pixel value 
is the median flux of the co-added ditherings normalized to the 
total exposure time and relative to the airmass and atmospheric 

transparency of the first frame in the dithering set. The final 

See: http://www.ls.eso.org/ | lasilla/sciops/2p2/E2p2M/WFI/CalPlan/fiat^kgd images were also trimmed to the overlapping area of 
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Table 1. Journal of the observations. 



Date 


Wrl riltei 


MJD 


Tcxp 


Air Mass 


Date 


Wri r liter 


MJD 


Tcxp 


Air Mass 


(d/m/y) 






(sec) 




(d/m/y) 






(sec) 




24/1 1/2001 


V/89 


52237.22128823 


3x10 


1.50 


04/03/2002 


V/89 


52337.20706880 


3x10 


1.22 




Rc/162 


52237.22634835 


3x5 


1.47 




Rc/162 


52337.21803382 


3x5 


1.25 




Ic/Iwp 


52237.23117713 


3x7 


1.44 




Rc/162 


52337.22950399 


3x25 


1.29 




Ic/Iwp 


52237.23529559 


3x25 


1.42 




V/89 


52337.23437902 


3x90 


1.31 




Rc/162 


52237.24074915 


3x25 


1.39 




V/89 


52337.24109513 


3x10 


1.34 




V/89 


52237.24882797 


3x90 


1.35 


05/03/2002 


V/89 


52338.21350087 


3x10 


1.24 




V/89 


52237.26303437 


3x10 


1.30 




Rc/162 


52338.21877655 


3x5 


1.26 




Rc/162 


52237.26868166 


3x5 


1.28 




Ic/Iwp 


52338.22356216 


3x7 


1.28 




Ic/Iwp 


52237.27821246 


3x7 


1.24 




Ic/Iwp 


52338.22763522 


3x25 


1.29 




Ic/Iwp 


52237.28227565 


3x25 


1.23 




V/S9 


52338.23909508 


3x90 


1.34 




Rc/162 


52237.28773647 


3x25 


1.21 


06/03/2002 


V/89 


52339.18778127 


3x10 


1.18 




V/89 


52237.29344333 


3x90 


1.20 




Rc/162 


52339.19297910 


3x5 


1.19 




V/89 


52237.30066930 


3x10 


1.18 




Ic/Iwp 


52339.20270214 


3x7 


1.22 




Rc/162 


52237.31212718 


3x5 


1.16 




Ic/Iwp 


52339.20714999 


3x25 


1.23 




Ic/Iwp 


52237.31721340 


3x7 


1.15 




Rc/162 


52339.21242581 


3x25 


1.25 


26/11/2001 


v/89 


52239.18065529 


3x10 


1.78 




V/89 


52339.21805674 


3x90 


1.27 




Rc/162 


52239.18562202 


3x5 


1.74 




V/89 


52339.22604426 


3x10 


1.30 




Ic/Iwp 


52239.19029892 


3x7 


1.69 




Rc/162 


52339.23099193 


3x5 


1.32 




Ic/Iwp 


52239.19447743 


3x25 


1.66 




Ic/Iwp 


52339.23578713 


3x7 


1.34 




Rc/162 


52239.19983479 


3x25 


1.62 




Ic/Iwp 


52339.23999762 


3x25 


1.35 




V/89 


52239.20535873 


3x90 


1.57 




Rc/162 


52339.24539390 


3x25 


1.38 




V/89 


52239.32447597 


3x10 


1.13 




V/89 


52339.25120870 


3x90 


1.41 




Rc/162 


52239.32957951 


3x5 


1.12 




V/S9 


52339.25913018 


3x10 


1.45 




Ic/Iwp 


52239.33441270 


3x7 


1.12 




Rc/162 


52339.26466470 


3x5 


1.49 




Ic/Iwp 


52239.33845778 


3x25 


1.11 




Ic/Iwp 


52339.26890995 


3x7 


1.51 




Rc/162 


52239.34398065 


3x25 


1.11 




Ic/Iwp 


52339.27256372 


3x25 


1.53 




V/89 


52239.35053232 


3x90 


1.10 




Rc/162 


52339.28297103 


3x25 


1.61 


23/01/2002 


V/89 


52297.19774832 


3x10 


1.10 




V/89 


52339.28856480 


3x90 


1.66 




Rc/162 


52297.20344695 


3x5 


1.10 




V/S9 


52339.29614328 


3x10 


1.73 




Ic/Iwp 


52297.20807586 


3x7 


1.09 




Rc/162 


52339.30145210 


3x5 


1.78 




Ic/Iwp 


52297.21229543 


3x25 


1.09 




Ic/Iwp 


52339.3060641 1 


3x7 


1.82 




Rc/162 


52297.21764270 


3x25 


1.09 




Rc/162 


52339.31534623 


3x25 


1.93 




V/89 


52297.22325254 


3x90 


1.09 


24/03/2002 


V/89 


52357.07361738 


3x10 


1.10 




V/89 


52297.23089695 


3x10 


1.09 




Rc/162 


52357.07891852 


3x5 


1.10 




Ic/Iwp 


52297.24111660 


3x7 


1.10 




Ic/Iwp 


52357.08412627 


3x7 


1.10 




Ic/Iwp 


52297.24521964 


3x25 


1.10 




IcAwp 


52357.08840132 


3x25 


1.10 




Rc/162 


52297.25143005 


3x25 


1.10 




Rc/162 


52357.09377871 


3x25 


l.ll 




V/89 


52297.25701984 


3x90 


1.11 




V/89 


52357.09933803 


3x90 


1.12 


11/02/2002 


V/89 


52316.14784279 


3x10 


1.10 




V/89 


52357.10984608 


3x10 


1.13 




Rc/162 


52316.15586182 


3x5 


1.09 




Rc/162 


52357.12092981 


3x5 


1.15 




IcAwp 


52316.16066498 


3x7 


1.09 




Ic/Iwp 


52357.12572673 


3x7 


1.16 




Ic/Iwp 


52316.16510996 


3x25 


1.09 




Ic/Iwp 


52357.13007714 


3x25 


1.16 




Rc/162 


52316.17086954 


3x25 


1.09 




Rc/162 


52357.13532178 


3x25 


1.17 




V/89 


52316.17698296 


3x90 


1.09 




V/89 


52357.16004606 


3x10 


1.24 




V/89 


52316.18480783 


3x10 


1.10 




Rc/162 


52357.16513120 


3x5 


1.26 




Rc/162 


52316.19440469 


3x5 


1.10 




IcAwp 


52357.16995618 


3x7 


1.27 




Ic/Iwp 


52316.19919322 


3x7 


1.10 




IcAwp 


52357.17390633 


3x25 


1.28 




Ic/Iwp 


52316.20328731 


3x25 


1.11 




Rc/162 


52357.18004707 


3x25 


1.31 




Rc/162 


52316.20891675 


3x25 


1.11 




V/89 


52357.18570452 


3x90 


1.33 




V/89 


52316.21545996 


3x90 


1.12 




V/89 


52357.99572072 


3x10 


1.13 


24/02/2002 


V/89 


52329.19058819 


3x10 


1.14 


25/03/2002 


Rc/162 


52358.00131936 


3x5 


1.13 




Rc/162 


52329.19609543 


3x5 


1.14 




Ic/Iwp 


52358.00615478 


3x7 


1.12 




Ic/Iwp 


52329.20090155 


3x7 


1.15 




Rc/162 


52358.01631727 


3x25 


1.11 




Rc/162 


52329.21181380 


3x25 


1.17 




V/89 


52358.02278763 


3x90 


1.10 




V/89 


52329.21777151 


3x90 


1.19 













the relative ditherings, so that only those sources falling in all 
ditherings are considered for the subsequent analysis. 

2.1.3. Photometric calibration 

Instrumental magnitudes were reported to the standard 
Johnson-Cousins system. To this aim, the Landolt's standard 



star field SA 92 (lLandoltll992h was observed in the VRI filters. 
By using the IRAF photcal package, we first performed the 
aperture photometry for the standard stars, obtaining their in- 
strumental magnitudes (vq, rg and /q) corrected for atmospheric 
extinction and normalized to the exposure time. Then, the 
transformation coefficients, namely zero point (ZP) and color 
term (c), from the WFI photometric system to the Johnson- 
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Table 2. Mean photometric calibration coefficients for our 
WFI observing run. 



Filter 




ZP 




c 


V 


0.130 


24.195+0.035 


0.071 


+0.037 


R 


0.096 


24.509±0.027 


0.078 


+0.031 


I 


0.082 


23.502+0.034 


0.325 


+0.044 



"'" AE = Atmospheric extinction coefficient. The mean values for La 
Silla have been adopted. 



Cousins standard system were determined by a linear fitting 
of the following equations: 



V ^Vi) + cv (vo - r()) + 2Pv 
Rc ^ ro + CR- (ro - /q) + ZPr 
Ic = io + c, ■ (ro - io) + ZP, 



(1) 
(2) 
(3) 



where V, Rc and Ic are the standard magnitudes of the 
Landolt's stars. The mean transformation coefficients deter- 
mined in our observing run are reported in Table |2]and are 
consistent with the mean values computed for La Sillo 

2.1 .4. The deep catalog extraction 

The sources extraction and photometry from each stacked im- 
age in each filter was pe rformed using the SExtractor tool by 
Bertin & Arnouts (1 1996 *). This data set will be used to study 
the variability of the IC 2391 low-mass members and will be 
presented in a future paper (Pagano et al., in preparation). The 
limiting magnitude achieved at the 3cr level using the single 
staked images is R ^20. At the age (~50 Myr) and distance 
(~150 pc) to IC 2391, this limit corresponds to very low-mass 
stars (MwO. 1 Mq) close to the hydrogen burning hmits or, per- 
haps, to massive BDs. 

One of the main goals of this work is the identification of 
new candidate substellar members of IC 2391, which are ex- 
pected to have R magnitude in the range 18.5-24 depending 
on their mass and the intervening reddening. Since our data 
set consists of repeated observations of the same field, we fur- 
ther combined the stacked images with the same exposure time 
and obtained, for each filter, a photometrically deeper frame 
with improved signal-to-noise ratio, suitable for our purposes 
(Fig. |2]l. The combination has been performed using the im- 
combine task under IRAF in the same configuration as ex- 
plained in Sect. 12.1.21 and the stacked images from our "long- 
time" series (see Sect. 12. II ). 17 images were combined in the V 
filter, 18 in the R filter and 15 in the / filter, for a total exposure 
time of 76.5 min, 22.5 min and 18.75 min, respectively. PSF- 
fitting photometry was performed on these deep er frames by 
using the IRAF daophot package (IStetsonlll987l) : this method 



allows a better deblending of sources with respect to aperture 
photometry and hence, the detection of the faintest neighbors to 
bright stars. Since we are interested in very low-luminosity ob- 




Fig. 2. Section of a single stacked image in the R-bsnd (upper 
panel) and the same region as it appears in the deeper frame 
obtained by combining 18 stacked images (lower panel). 



jects, we set a very low detection threshold (cr=3) and then se- 
lected point-like objects from extended sources (e.g. galaxies), 
saturated objects, and other spurious detections by using the 
daophot morphological parameters. In Fig.|3]the internal pho- 
tometric errors of all the detected point-like sources are plotted 
against magnitude for aU the available filters; the relative expo- 
nential fits are over-plotted. Table |3] reports the saturation lim- 
its and the limiting magnitudes achieved at the 3cr level from 
these deeper frames. These limits are suflicient for the detec- 
tion of young stars and BDs with mass from ~1 down to 
-0.03 Mq at the age and distance of IC 2391. 

Our optical catalog contains 60596 point-like sources with 
VRI magnitudes. A merged catalog with position and photom- 
etry in each optical filter, as well as in th e near-IR {JHKs 
magnitudes from 2MASS; Cutri et al. 2003b . was finally pro- 
duced and used for the selection of cluster members (Sect.O. 
A matching radius of l"was used to merge the WFI and the 
2MASS catalogs; this value was set by taking the astrometric 
accuracy of both catalogs into account. 



See http://www.ls.eso.org/lasilla/Telescopes/2p2T/E2p2M/WFI/zero 



2.1 .5. Completeness 

The completeness of our "deep" catalog was derived empiri- 
cally. We used the IRAF addstar package to add test stars into 
the frames and then the same PSF-fitting procedure as for the 
real objects (Sect. in~4l i to determine what fraction of these ar- 
tificial stars were recoverable as a function of magnitude. This 
fraction provides a measure of the completeness. We inserted 
500 artificial sources in each of the three deep frames in the 
VRI filters; this number is statistically significant for our pur- 
pose and, at the same time, does not alter the crowding in the 
images. The profile for the artificial sources was generated by 
iptamg'.the same PSF model as used for the source extraction; 
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Fig. 3. Photometric errors versus magnitudes and relative ex- 
ponential fits for all the point-like sources detected in the 
30x30 square arcmin surveyed area in IC 2391. Magnitudes 
and errors are from the deep frames obtained as discussed in 
Sect. l2T4l 

Table 3. Saturation limits and limiting magnitudes at 3cr level 
(see Sect. l2.1.4] l. In the fourth column the corresponding limit- 
ing mass for cluster members (i.e. assuming an age between 30 
and 50 Myr and a distance of 150 pc) derived from the evo- 
lutiona ry models by iBaraffe et al.l (Il998i) and Ichabrier et aL 
(12000 ') is reported; in the last column the magnitude cor- 
responding to the Hydrogen burning limit (~0.07 M©) for 
IC 2391 members is also given for comparison. 



Filter 


Mag Sat. 


Mag 3cr 


Mass 3cr 
(Mo) 


PHB' 


V 


10.5 


24.7 


0.03 


20.0 


Rc 


10.0 


23.7 


0.03 


18.5 


Ic 


9.5 


23.0 


0.02 


17.0 



t PHB = Photometric Hydrogen Burning limit for IC 2391 members. 



the positions of the artificial objects were randomly distributed 
over the entire area of the mosaic, and their magnitudes range 
uniformly between the detection and the saturation limits in 
each band (Table [3]). Fig. |4] shows the fraction of recovered 
artificial objects, i.e. the completeness (C), as a function of 
magnitude for each filter We estimate a 100% completeness 
level down to V ~20, R ~19.5 and / «I8.5; however, objects 
down to our detection limits at 3cr level (Table I2.1.4I I are re- 
covered with a completeness level better than ~8 0%. Using the 



theore tical isochrones ed evoluti onary tracks by IBaraffe et al 
(Il998h and lChabrier et all (l2000l) for the WFI-Cousins system 
(see Sect. 3.1 by ISpezzi et al.ll2007h and considering the age 
and distance of IC 2391, we estimate that the 100% and 80% 




Mognitude 

Fig. 4. Completeness (C) plot for extraction of artificial stars 
from our "deep" mosaics (see Sect. 12.1.41) for the three optical 
bands used in this work. 



completeness levels correspond to objects with Ma!0.05 Mq 
and Mx;0.03 M©, respectively. Interstellar extinction effects are 
expected to be negligible (see Sect.[T]i. 

2.2. Spitzer observations and data reduction 

IRAC dFazio et all l2004 observations of IC 2391 were ob- 



tained on 30 December 2006. The 12s high-dynamic-range 
mode was used to obtain two frames in each position, one with 
0.4 sec exposure time and one with 10.4 sec. The observation 
of each field was repeated three times with a small offset, pro- 
viding 31.2 sec integration time for each position. The frames 
were processed using the SSC IRAC Pipeline vl5.05, and mo- 
saics were created from the basic c ahbrated data (BCD) fra mes 
using a custom IDL program (see Guterm uth et al. I (Egos') for 
details). Aperture photometry on these images was carried out 
using Phot Vis version 1.10, w hich is an IDL-GUI bas ed pho- 
tometry visualization tool. See lGutermuth et alJ (120041) for fur- 
ther details on PhotVis. The radii of the source aperture, and of 
the inner and outer boundaries of the sky annulus were 2.4, 2.4 
and 7.2 arcsec, respectively. The calibration used large aperture 
measurements of standard stars. The zero point magnitudes of 
the calibration were 19.6642, 18.9276, 16.8468, and 17.3909 
for channel 1, 2, 3, and 4, respectively. Aperture corrections 
of 0.21, 0.23, 0.35 and 0.5 mag were applied for channels 1, 
2, 3, and 4 to account for the differences between the aperture 
sizes used for the standard stars and for the IC 2391 photom- 
etry. The sensitivity of the IRAC observations is sufficient for 
the detection of objects with Ks <14.5, which corresponds to 
the hydrogen-burning limit for cluster members. 

MIPS observations in IC 2391 and the relative data 
reduction and cali bration have been already described by 
SiegIeretalJ(l2007h . 

The standard fluxes at each Spitzer pass-band were derived 
from the observed magnitudes using the absolute flux call- 
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bration constants bv lReach et alj ( l2005b and lEngelbracht et al 
( l2007h for IRAC and MIPS/24jum bands, respectively. 



3. Selection of new cluster members 

The extraction of member candidates in open clusters generally 
consists in selecting objects whose position in color-magnitude 
diagrams (CMDs) is located abov e the ZAMS shifted to the 
distance of the cluster (see, e.g., iBarrado v Navascues et al. 
2001 at lLeggettl ll992t lLodieull2005b . The selection of IC 2391 
member candidates was carried out using a similar criterion, re- 
fined by exploiting estimates of the cluster proper motion and 
the photometry of previously known cluster members available 
from the literature. Our selection procedure includes the three 
steps which are now described. 



3.1. Determination of stellar parameters 

We first assumed that all the point-like objects detected in 
our images above 3cr level are all at the distance of IC 2391 
and simultaneously estimated their visual extinction (Ay), ef- 
fective temperature (Teff) and stellar radius (R*) by fitting 
a grid of stellar photosphere models to their observed spec- 
tral energy distribution (SED). The observed SEDs were de- 
rived by merging the VRI WFl-Cousins photometry with the 
2MASS JHKs photometry; the standard fluxes at each pass- 
band were derived from the observed magnitud e s usin g the 
absolute flux calibration constants by 'Johnson' ("1965') and 
Cousins d 19761) for the optical bands and Cutri et al. (2003) for 
the 2MASS bands, respectively. The grid of reference SEDs 
was const ructed by using t he Ne xtG en and Stardusty rn odel 
spectra bv lHauschildt et al.l (Il999b and lAllard et al.l (|2000|) and 
the WFI and 2MASS fi lter transmission curves, as prescribed 
by ISpezzi et al.l OOOTh . Our grid spans the Teff range 1700- 
10000 K and assumes a surface gravity \ogg = 4.5, as ex- 
pected for low-mass objects at the age of IC 2391. The refer- 
ence SEDs were scaled to the distance to IC 2391 (i.e. 150 pc) 
and fitted to t he observed SED o f each object following the 
prescription by Spezzi et al. (2007) and adopting the extinction 
law by CardelUet al.. (.1989^1 The Ay, Te« and R* values of the 
reference SED minimizing the of the fit correspond to the 
best approximation of the actual stellar parameters. The stellar 
luminosity (L*) is derived as = AncrRi,2T'^ff, where cr is 
the Stefan-Boltzmann constant. For more details on the SED 



fitting procedure we defer the reader to ISpezzi et al.l (l2007h . 
The goodness of our fits can be appreciated from the examples 
shown in Fig.|8] 



As demonstrated by ISpezzi et al.l (12007 ). this parametriza 



tion method turns out to be accurate within about 250 K in Tgff 
and 1.5 mag in Ay for objects with spectral type later than K5 
(i.e Tj,// <4500 K). Such uncertainties translate into an av- 
erage uncertainty of 0.15 dex in the logarithm of luminosity. 



^ We also ch ecked that using other extinction laws, such as those 
by ISavage & M athis (1979) and Weingartner & Draine (2001J does 
not affect significantly our results. We adopt the lCardelli et all j 19891) 
extinction law because its analytic form is more functional for our 
purposes. 



This accuracy in the derived physical parameters is sufficient 
for candidate selection purposes. Moreover, on the basis of the 
saturation limits of our survey (Table |3j, we expect to recover 
only the low-mass cluster members (Rc >10 mag, i.e. M<1 Mq 
at the age and d istance of IC 2391), for which the parametriza- 
tion method by Spezzi et al.l (12007!) can be confidently applied. 
Note also that possible IR excesses at 2MASS wavelengths due 
to the presence of hot circumstellar dust does not appear to be 
a major issue in our fits. Indeed, the age of IC 2391 is consis- 
tent with both theoretical and observational tirn^^cales of disk 
dissipation and terrestrial planet for mation dChamberslboOll: 
Kleine et al.l2002l:ISiegler et al.l2007b : di us, following the clas - 



sification by Lada| (jl987|) as extended by Greene et al. ( 1994 ) 



the cluster population is expected to consist mainly of class III 
objects with no prominent IR excess at near-IR wavelengths. 
Moreover, the same parametrization procedure was applied to 
40 confirmed members of IC 2391 whose spectral type (later 



the literature (IBarrado v Navascues et al. 


1999i boOlb. 2004; 


Cutri et al.! 20031 IPatten & Simon 1996 




Patten <& PavlovskJ 


1999;,Rolleston et al.a 997;, Simon & Pattenia998'). Our fitting 



procedure yields for these objects Teff and Ay values consis- 
tent with those derived from spectroscopy within -250 K and 
-0.7 mag, respectively. This further supports the reliability of 
our parametrization method. 

3.2. Selection of cluster member candidate from the 
HR diagram 

Having a reliable estimate of Teff and for all the sources 
in our catalog with Teff <4500 K (i.e. 2395 sources), we 
constructed the Hertzsprung-Russell (HR) diagram for these 
sources and for the 40 confirmed members of IC 2391 men- 
tioned above (Sect. 13.1b . Mass and age of these objects can 
be now determined by comparison of their location on the HR 
diagram with a suita ble set of evo l utiona ry tr acks, which in 



our case are those by Baraff e et al.l (119981) and IChabrier et al 



(2000). We find that the isochronal age of the 40 low-mass clus- 
ter members ranges between 10 and 80 Myr, with a mean value 
of 20-30 Myr (Fig. |5]i; this isochronal age is in fair agreement 
with the cluster age estimated on the basis of the lithium de- 
pletion te chnique and the main-sequence is o chron e fitting (30- 
50 Myr; 'Barrado v Navascues et al.' 'l999', 12004 iBasri et al 



[1996; Mermilliod 1981; Stauffer et al. 1998). However, we ob- 
serve an age spread between a few Myr and 100 Myr (see also 
TablelU which is larger than that reported in the literature. This 
larger age spread mainly arises from our indetermination on ef- 
fective temperature and luminosity (Sect. 13.11 1; indeed, stellar 
parameters inferred from photometry alone are far more inac- 
curate than those derived from spectroscopy, though still useful 
for candidate selection purposes. 

New low-mass cluster members are expected to be found in 
the same locus of the HR diagram. Thus, we performed a first- 
level selection of member candidates by picking up those ob- 
jects whose isochronal age is <100 Myr; this selection yielded 
90 candidates. Clearly, we also find a large number of source in 
our catalog for which we cannot assign an age because their lo- 
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Fig. 5. HR diagram for 40 low-mass confirmed members of 
IC 2391 (stars), the sample of optical sources whose isochronal 
age is <100 Myr (open squares; see Sect. 13.21 1 and the final 
sample of 29 member candidates selected in this work (filled 
squares). The continuous and dashed lines are the evoluti onary 



track s and isochrone by Baraffe et al. ( 1998) and Chabri er et al 



(12000 ). The average error bars on effective temperature and lu- 
minosity are shown in the top-right. 



cation on the HR diagram is inconsistent with the pre-main se- 
quence (PMS) locus when the distance of IC 2391 is assumed; 
these sources are field objects unrelated to the cluster. 

The advantage of using this approach, as compared to fre- 
quently used membership criteria based on CMDs, lies in the 
simultaneous use of the photometric information from the V to 
the Ks bands; indeed, this wavelength range samples the peak 
of the photospheric emission of cool, low-mass objects, and en- 
ables a highly reliable determination of the stellar parameters 
(Comeron et al., in preparation). 

3.3. Refining the sample using proper motion 
measurements 

Fig. |6] show the Rc vs. ( V - Rc) and Ic vs. (R - I)c CMDs 
for all the point-like optical sources detected in IC 2391 above 
3(T level and the 90 candidates selected from the inspection 
of the HR diagram. The lines in these plots are the 1 and 



50 Myr theoretical NextGen isochrones bv lBarafFe et al.l ( 119981) 
for effective te mperatures high er than 2500 K and the DUSTY 



isochrones by IChabrier et al. I (EOOQ) for lower temperatures; 
these isochrones were transformed i nto th e WFI-Cousins sys- 
tem, as described by ISpezzi et al.l (120071) . and scaled to the 
IC 2391 distance of 150 pc. No correction was made for the 
mean reddening towar ds IC 2391, which is estimated close 
to zero (E(B-V)^0.006; |Patten & Simonll 19961) . Both diagrams 
show a wide sequence of candidates with 12 < Rc <15 mag 
and 11< Ic 5 14 mag, which is likely due to the merging of 
the cluster and field star sequences. Indeed, since the galactic 
latitude of IC 2391 is about -7°, the sample of photometrically 
selected candidates is inevitably contaminated by field stars. 



I I I 1 1 I I 1 1 I 1 1 I I 1 1 I I 1 1 I 1 1 I I 1 1 I 1 1 
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Fig. 6. Rc vs. (V - Rc) and Ic vs. (R - I)c CMDs for the 
point-like sources in our WFI survey (small gray dots). The 
lines represent the 1 and 50 Myr theoretical isochrones derived 
as explained in Sect. 13. 3| shifted to the distance modulus of 



IC 2391 (5.88 mag; iForbes et al.l 1200 Ih . The open squares 



represent the sample of optical sources whose isochronal 
age is <100 Myr (Sect. I3.2l i. while the filled squares mark 
the final sample of 29 member candidates (Sect. I3.3I I. The 
stars are the confirmed cluster m embers, whose photometry 



was re tr ieved from the l it eratur e dBarrado v Navascues et al 



RoUeston et al.l 



1997t 



2001bl; IPatten & SimonI Il996l; IPatten & Pavlovskil Il999 



Simon & Patten 1998h. Error bars are 



also drawn; where not visible, they are smaller than the symbol 



In order to weed out stellar contaminants, proper mo 
tions of cluster member candidates are generall y computed and 



Casewell et al 



compared to those of field objects (see, e.g., 
,2007,) . We applied this proper motion analysis to the candi- 
dates in our sample using two proper motion catalogs: the all- 
sky Naval Observator y Merged Astrometric Dataset (NOMAD 



Zacharias et al ]|2005h . reporting proper motions for stars down 



to ~16 mag with an average accuracy of 6 mas/yr, and 
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Table 4. Optical and near-IR photometry for the IC 239 1 members and candidates of IC 239 1 discussed in this work. The objects 
with identification number (ID) from 1 to 29 are the new candidates selected in this work (Sect.|3li; three of these candidates show 
IR excess emission at 24//m and are marked by an asterisk. The objects with ID from 30 to 37 are the additional 8 cluster members 
and candidates showing IR excess emission discussed in Sect. |6] 



ID 


Designation 


Status^ 


Ref." 


„ . t 
'■'■'^■72000 
{hh:mm:ss) 


(dd:mm:ss) 




V" 




,t 








1 


WFI J08403964-5250429* 


HR-C 




08:40:39.643 


-52:50:42.86 


- 


21.24±0.09 


19.82+0.05 


17.83+0.07 


15.54+0.05 


14.86+0.08 


14.52+0.09 


2 


WFI J08404512-524852I 


HR-C,PM-C 




08:40:45.118 


-52:48:52.09 


17.10 


16.60±0.07 


15.39+0.03 


13.97+0.04 


12.23+0.02 


11.65+0.02 


11.42+0.02 


3 


WFI J08404578-5301052 


HR-C 




08:40:45.780 


-53:01:05.23 


- 


14.50+0.07 


13.36+0.03 


12.63+0.04 


11.24+0.03 


10.43+0.04 


10.19+0.05 


4 


WFI J084055 18-5257227 


HR-C 




08:40:55.176 


-52:57:22.68 


15.52 


14.59±0.07 


13.39+0.03 


12.58+0.04 


10.94+0.03 


10.10+0.03 


9.84+0.02 


5 


WFI 108410806-5300053°'' 


HR-C,X 


a 


08:41:08.057 


-53:00:05.26 


- 


18.02+0.08 


16.82+0.06 


15.27+0.06 


13.28+0.03 


12.67+0.06 


12.45+0.05 


6 


WFIJ084I08 14-5255027 


HR-C.PM-C 




08:41:08.141 


-52:55:02.71 


14.30 


14.69+0.07 


13.45+0.03 


12.64+0.04 


10.98+0.04 


10.08+0.04 


9.81+0.03 


7 


WFI J084I0886-5248387* 


HR-C.X 




08:41:08.856 


-52:48:38.74 


18.77 


18.74+0.07 


17.46+0.04 


15.66+0.06 


13.52+0.03 


12.90+0.02 


12.65+0.03 


8 


WFIJ084I0937-5302136 


HR-C,X 




08:41:09.374 


-53:02:13.60 


17.51 


17.58+0.07 


16.38+0.04 


14.74+0.05 


12.73+0.03 


12.11+0.03 


11.91+0.03 


9 


WFI J0841 1 865-5258549™ 


HR-C,X 


a 


08:41:18.650 


-52:58:54.88 


19.61 


18.83+0.07 


17.64+0.03 


16.12+0.05 


14.14+0.03 


13.58+0.03 


13.23+0.04 


10 


WFIJ084I2229-530447I 


HR-C,X 




08:41:22.291 


-53:04:47.14 


17.77 


17.16+0.07 


15.98+0.03 


14.43+0.05 


12.47+0.02 


11.85+0.02 


11.58+0.02 


11 


WFIJ08413933-5252565 


HR-C 




08:41:39.326 


-52:52:56.53 


15.76 


14.93+0.07 


13.81+0.03 


13.09+0.04 


11.60+0.02 


10.85+0.03 


10.65+0.02 


12 


WFIJ08415360-5257558 


HR-C,X 




08:41:53.602 


-52:57:55.80 


18.12 


17.44+0.07 


16.24+0.03 


14.72+0.05 


12.79+0.03 


12.17+0.02 


11.93+0.03 


13 


WFI J08422460-5259575 


HR-C.PM-C 




08:42:24.605 


-52:59:5748 


13.71 


12.71+0.05 


11.91+0.04 


11.48+0.05 


10.54+0.03 


10.05+0.03 


9.92+0.02 


14 


WFII0842307I-5257348 


HR-C,PM-C,X 




08:42:30.715 


-52:57:34.78 


14.95 


14.63+0.07 


13.49+0.03 


12.73+0.04 


11.39+0.02 


10.71+0.02 


10.52+0.02 


15 


WFI J0842473I-5309235 


HR-C.X 




08:42:47.314 


-53:09:23.54 


17.80 


20.36+0.08 


18.91+0.05 


16.79+0.07 


14.20+0.05 


13.64+0.05 


13.24+0.04 


16 


WFI J08424850-5303I87 


HR-CX 




08:42:48.502 


-53:03:18.72 


19.07 


18.55+0.07 


17.26+0.04 


15.63+0.05 


13.59+0.02 


12.99+0.02 


12.68+0.03 


17 


WFI J08425873-5305546 


HR-C,PM-C,X 




08:42:58.733 


-53:05:54.64 


17.56 


16.55+0.07 


15.33+0.03 


13.87+0.04 


12.01+0.02 


11.38+0.02 


11.13+0.02 


18 


WFI 108425988-5305074 


HR-C 




08:42:59.880 


-53:05:07.44 




20.65+0.08 


19.19+0.05 


17.17+0.07 


14.77+0.04 


14.16+0.04 


13.81+0.06 


19 


WFI J08430294-5302259 


HR-C 




08:43:02.938 


-53:02:25.91 


15.07 


13.97+0.07 


12.87+0.03 


12.10+0.04 


10.60+0.02 


9.79+0.02 


9.59+0.02 


20 


WFIJ08431228-525141I 


HR-C 




08:43:12.283 


-52:51:41.11 




22.85+0.10 


21.32+0.07 


18.92+0.09 


16.18+0.09 


15.53+0.13 


15.45™ 


21 


WFIJ0843 1 860-5306241 


HR-C,PM-C 




08:43:18.602 


-53:06:24.12 


13.97 


1341+0.06 


12.43+0.03 


11.80+0.04 


10.52+0.02 


9.83+0.02 


9.67+0.02 


22 


WFI J08433366-5259549 


HR-C 




08:43:33.658 


-52:59:54.89 




22.71+0.12 


20.88+0.07 


18.55+0.08 


15.80+0.08 


15.05+0.09 


14.75+0.11 


23 


WFI J08433843-5246379 


HR-C 




08:43:38.431 


-52:46:37.92 


18.09 


16.96+0.08 


15.49+0.03 


14.36+0.04 


12.38+0.02 


11.36+0.02 


11.06+0.02 


24 


WFI J08433939-5307563™ 


HR-C 


b 


08:43:39.389 


-53:07:56.32 


18.38 


18.08+0.07 


16.87+0.04 


15.20+0.05 


13.14+0.02 


12.54+0.03 


12.28+0.02 


25 


WFI J08433945-5248576 


HR-C 




08:43:39.451 


-52:48:57.56 


19.74 


18.79+0.07 


17.60+0.03 


16.16+0.04 


14.38+0.03 


13.76+0.03 


13.46+0.05 


26 


WFI J08434640-5257036 


HR-C.PM-C 




08:43:46.399 


-52:57:03.56 


14.48 


13.73+0.06 


12.77+0.03 


12.12+0.04 


10.86+0.03 


10.23™ 


10.01™ 


27 


WFII0843475I-5250I36 


HR-C.PM-C 




08:43:47.508 


-52:50:13.63 


16.01 


15.92+0.07 


14.59+0.03 


13.64+0.03 


11.88+0.02 


10.91+0.02 


10.65+0.02 


28 


WFI J08434909-5307269 


HR-C 




08:43:49.090 


-53:07:26.94 


19.51 


18.64+0.07 


17.48+0.03 


16.00+0.04 


14.13+0.03 


13.55+0.03 


13.28+0.04 


29 


WFI J08440199-5310356 


HR-C.PM-C 




08:44:01.992 


-53:10:35.62 


15.20 


14.95+0.06 


13.94+0.03 


13.24+0.04 


11.65+0.04 


11.09+0.05 


10.69+0.03 


30 


WFI I08404683-5302487*™ 


M 




08:40:46.826 


-53:02:48.73 




19.25+0.06 


18.36+0.03 


17.80+0.04 


16.53+0.21 


15.81+0.28 


15.95"'- 


31 


WFI J084I0019-5249088* 


IRC 




08:41:00.190 


-52:49:08.76 


17.53 


16.71+0.07 


15.60+0.03 


14.92+0.04 


13.51+0.03 


12.77+0.03 


12.56+0.03 


32 


WFI J08412176-5247289»''^ 


M 


c 


08:41:21.756 


-52:47:28.90 


20.10 


19.33+0.06 


18.24+0.03 


17.65+0.04 


16.43+0.11 


15.60+0.14 


15.25™ 


33 


WFI J08412550-5300154* 


IRC 




08:41:25.502 


-53:00:15.44 




18.73+0.05 


17.92+0.04 


17.49+0.05 


16.35+0.11 




15.60™ 


34 


WFI J08413601-5309274*™ 


M 


d,c 


08:41:36.010 


-53:09:27.40 


19.37 


18.52+0.07 


17.27+0.04 


15.56+0.05 


13.43+0.02 


12.89+0.03 


12.61+0.03 


35 


WFI J08413888-5259492* 


IRC 




08:41:38.885 


-52:59:49.24 


17.22 


17.20+0.06 


16.30+0.03 


15.73+0.04 


14.56+0.02 


14.16+0.03 


13.95+0.06 


36 


WFI J08414904-5255192* 


IRC 




08:41:49.039 


-52:55:19.20 


19.56 


18.52+0.06 


17.46+0.03 


16.77+0.04 


15.19+0.05 


14.61+0.06 


14.47+0.10 


37 


WFI 108415780-5252140*''^ 


M 


f 


08:41:57.799 


-52:52:13.98 


14.74 


13.76+0.06 


12.74+0.03 


12.16+0.04 


11.09+0.02 


10.40+0.02 


10.27+0.02 



M = confirmed cluster mem ber; HR-C = memb er candidate according to the location on the HR diagram; PM-C = proper motion member 
can didate; X = X-ra y sourc es jMarino et al.ll2005h.; IR-C: you ng stellar obje ct ca ndidate accor d ing to the IR properties (Sect. 6.1b . 

* a) iRolleston et al] h997h ; b) IPatten & Pavlovskil h999l) ; c) lOodd ( l2004h ; e) | Marino et al ] bOOSh ; f) ISiegler et al.l ( 1200717 and references 
therein. 

^ From the WFI survey presen ted in this work. 

" From the NOMAD catalog (Zacharias e t alj2005h . Photometric uncertainties are of the order of 0.2-0.3 mag. 

* From the 2MASS catalog (Cutri et al. 200^. 

* Object showing IR excess emission at 24pm. 

Flux upper limit. 

''^ Other Designation. ID 5: RB729, ID 9: RB955, 2MASSI0841 1862-5258546. ID 24: PP9, 2MASSI08433935-5307556. ID 30: D97, 
2MASS J08404679-5302481. ID 32: D112, 2MASSJ08412176-5247291. ID 34: CTIOlOO, 2MASSJ08413595-5309268. ID 37: V368Vel, 
SMY29, PMM4636, GSC08569-01100, VXRPSPC41, SHJM9, VXRHRI19, #29. 



the catalog recently presented bv lPlatais et alj (120071) . who ob- 
tained proper motions for 6991 stars down to V«il3-16 over 
a ~9 deg^ area around IC 2391 with a precision of about 
1.7 mas/yr. These two catalogs report consistent proper mo- 
tions for the confirmed cluster members, however the NOMAD 
catalog is photometrically deeper and hence, we found a much 
larger number of matches with our candidate list with respect 
to the Platais et al. 1 (l2007h catalog Thus, we used the proper 
motions from NOMAD, where 70 candidates out of 90 in our 



original list (Sect. [J!2l ) have proper motion estimates. Because 
of the low precision, we cannot calculate a proba bility of mem 



bership by the standard method as described in Mor aux et al 



(2003). Consequently, we only selected objects to be proper 
motion members if their proper motions fall within 3cr of 
the cluster mean proper motion ipmiRA) x-24 mas/yr and 
pm(DEC) ~+24 mas/yr; Dodd 2004; Kharchenko 2001). This 
criterion yielded only 9 proper motion member candidates out 
of 70. Adding these objects to the 20 fainter candidates whose 
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proper motion are not available, we end up with a sample of 29 
new member candidates of IC 2391. Of course, we still expect 
a certain level of contamination, especially among the candi- 
dates lacking proper motion estimates; this issue is discussed 
in Sect.m 



3.4. Statement on the selection results 

Our selection method allowed us to identify 29 new candidate 
members of IC 2391: 9 selected by photometry and astrom- 
etry, and 20 selected by photometry, for which no PM data is 
available. In Table|4]and Table|5]we report the optical and near- 
IR photometry for these candidates and their stellar parameters 
as explained in Sect. 13.11 the object designat ion follows the 
scheme already adopted by Spezzi et al. ( 20071) . A few of these 
objects coincide with member candidates selected in pre vious 
photo metric surveys (Patten & Pavlovski 1999; Rolleston et al.l 
1997h . as indicated in Table |4] (column 5). Moreover, 10 of 
our candidates, i.e. 34%, are X-ray sources recently detected 
with XMM-EPIC by Marino et al. (2005); these 10 objects are 
marked with an "X" in Table|4](column 4). Note that the XMM- 
EPIC pointing in IC 2391 is centered at R.A.k 8''42'«00' and 
Decs; -53''00'36" and the instrument field of view has a di- 
ameter of 30 arcmin. Thus, t he great maiority of o ur candidates 
falls in the area observed by iMarino et alj ( l2005h and the rea- 
son why the 66% of them were not detected must be mainly 
ascribed to their faint X-ray emission (i.e. below the XMM- 
EPIC sensitivity limits: lEhle etalJl2008l) . 

In addition, we correctly recovered 32 previously known 
cluster members. We note that these 32 objects represent the 
~40% of the IC 2391 known population falling in the area ob- 
served with WFI; the remaining 60% are either saturated in our 
WFI images or do not have complete VRcIcJHKs photomet- 
ric data-set in our merged catalog, so we could not apply our 
selection criteria. However, using t heir ma gnitudes and proper 
motions from the t he literature fSarrad o v Navascues et al. 
1999', '2001b', '2004'; 'Cutri et al.' 120031: Patten & Simon 199& 
Patten & Pavlovski 1999; Platais et al.ll2007; Roll eston et all 
1997t lSimon & Patten 1998; Zacharias et al. 2005), they would 
fall within our selection (see Fig.|5]and Fig.|6]l. 

All these arguments support the reliability of our selection 
method. 



3.5. Multiple visual systems 

In order to identify possible multiple systems, we have visually 
inspected the 32 IC 2391 members falling in the area observed 
with WFI and the proposed 29 candida tes in all the available 
image s. Following the prescription by iReipurth & Zinnecker 
(119931) . we perform this search for visual neighbors in the sep- 
aration range 1.2-5 arcsec. The lower limit of this range was 
fixed considering the typical seeing of our images (Sect. 12. It . 
The upper limit is imposed to avoid contamination from fore- 
and background stars. This cut-off is somewhat arbitrary; we 
may justify this choice considering that the typical star-to-star 
separation on our images is ~5 arcsec and, when a visual com- 



panion is found within this distance, it mostly lies conspicu- 
ously close to the primary in an otherwise rather empty space. 

22 objects in our sample are found to have neighbors within 
less than 5 arcsec. However, a nearby star is accepted as a pos- 
sible companion only if it shows indications of youth. 2 of 
these neighbors have complete VRcIcJHKs photometric data- 
set but do not fulfill our membership criteria (Sect.O and hence 
they were rejected. The remaining 20 neighbors are too faint 
and do not have near-IR detection in 2MASS, so the only pos- 
sible indication of youth comes from their position on the op- 
tical CMDs. We found that 7 out of these 20 neighbors have 
optical magnitudes and colors consistent or marginally consis- 
tent with those of IC 2391 member^. These 7 objects, together 
with their possible primary companion and the relative angular 
separation, are reported in Table |6l Note that at the distance 
of IC 2391 (~150 pc), the range of angular separation cov- 
ered by these possible visual binaries corresponds to a linear 
separation of 350-700 AU, which is reasonable for physically 
bound young low-mass stars (Reipurth & Zinnecker 1993). For 
the member candidate ID 6 two neighbors were found within 
the 5 arcsec radius; both these neighbors have optical magni- 
tudes and colors consistent with membership to the cluster and 
hence this might be a triple system. 

4. Estimating the contamination 



In Sect. 13.31 we adopted the proper motion membership cri- 
terion to weed out stellar contaminants in the magnitude 
range Rc <16 mag. We could not perform this analysis for 
fainter candidates because their proper motion are unknown 
and hence, we still expect a certain level of contamination by 
field stars in our candidates sample. 

To estimate the degree of contamination we have used op- 
tical imaging of a nearby field available from the ESO archive. 
This field was observed with WFI in the VRI filters and is cen- 
tered at R.A.=08'':43'":53'' and Dec.=-52":47"':33', i.e. about 
0.9 deg north-east from the center of IC 239 1; since the cluster 



has an apparent diameter of about 50 arcmin (iKharchenko et al 



12005'), we expect this area to be more or less devoid of cluster 
members and we can use it as a control field. The exposure 
time of the control field images is 600 sec in all the filters; this 
implies that all the source brighter than R ^15 are saturated, 
while sources with 15< R <21 are detected with a sigma-to- 
noise ratio better than 10. Thus, using these data, we can es- 
timated the level of contamination at least in the photometric 
range spanned by our candidates with no proper motion deter- 
mination. By using the same procedure outlined in Sect. [3] we 
found 12 objects in our control field satisfying our photometric 
cluster membership criteria. Therefore, we estimate that the 17 
objects with 15< R <21 cataloged by us as candidate cluster 
members (TablelH still have a probability up to ~70% of being 
spurious. 

In order to further verify the amount of expected contam- 
ination in all the photometric range covered by our candi- 
dates, i.e. 12< Rc <21, we have also used a different ap- 
proach. Nearly all members of IC 2391 are located in the re- 



This membership criterion is explained in Sect. 13.3 
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Table 5. Estimates of the stellar parameters for the IC 2391 members and candidates of IC 2391 discussed in this work. The 
objects with identification number (ID) from 1 to 29 are the new candidates selected in this work (Sect. [3); three of these 
candidates show IR excess emission at 24jum and are marked by an asterisk. The objects with ID from 30 to 37 are the additional 
8 cluster members and candidates showing IR excess emission discussed in Sect. |6] 



ID 


Designation 










Mass* 


Age* 






(mag) 


(K) 


(Lo) 


(Ro) 


(Mo) 


(Myr) 


1 


WFI J08403964-5250429* 


1.63±0.13 


3000 


0.003±0.001 


0.19±0.04 


0.08±0.02 


56±20 


2 


WFIJ084045 12-5248521 


0.01 ±0.06 


3400 


0.038±0.017 


0.56±0.13 


0.25±0.05 


14±5 


3 


WFI J08404578-5301052 


0.97±0.07 


4200 


0.164±0.075 


0.76±0.18 


0.75±0.10 


45±20 


4 


WFI J084055 18-5257227 


1.80±0.06 


4200 


0.256±0.118 


0.86±0.22 


0.85±0.10 


32±10 


5 


WFI J08410806-5300053 


0.19±0.10 


3100 


0.015±0.007 


0.43±0.10 


0.07±0.02 


6±3 


6 


WF1J084108 14-5255027 


1.96±0.07 


4200 


0.261 ±0.1 20 


0.97±0.22 


0.90±0.10 


25±8 


7 


WFI J08410886-5248387* 


0.50±0.08 


3000 


0.012±0.006 


0.42±0.09 


0.07±0.02 


6±3 


8 


WFIJ08410937-5302136 


0.26±0.07 


3100 


0.024±0.01 1 


0.54±0.12 


0.13±0.06 


7±2 


9 


WFIJ0841 1865-5258549 


0.05±0.08 


3100 


0.007±0.003 


0.28±0.06 


0.11 ±0.04 


32±10 


10 


WFIJ08412229-5304471 


0.03±0.07 


3100 


0.031±0.014 


0.61±0.14 


0.11 ±0.04 


4±2 


11 


WFIJ08413933-5252565 


0.42±0.06 


4000 


0.093±0.043 


0.64±0.15 


0.62±0.08 


91 ±40 


12 


WFIJ08415360-5257558 


0.21±0.07 


3200 


0.023±0.011 


0.50±0.12 


0.20±0.04 


18±6 


13 


WFI J08422460-5259575 


0.07±0.06 


4800 


0.28B±0.133 


0.78±0.18 


0.85±0.10 


74±30 


14 


WFIJ08423071-5257348 


0.05±0.06 


4000 


0.102±0.047 


0.67±0.15 


0.70±0.10 


45±20 


15 


WFIJ0842473 1-5309235 


1.74±0.10 


2900 


0.008±0.004 


0.35±0.08 


0.05±0.02 


6±3 


16 


WFI J08424850-5303187 


0.03±0.07 


3000 


0.011±0.005 


0.38±0.09 


0.08±0.02 


9±4 


17 


WFI J08425873-5305546 


0.22±0.06 


3300 


0.049±0.023 


0.68±0.16 


0.25±0.05 


9±4 


18 


WFI J08425988-5305074 


1.74±0.10 


3000 


0.005±0.002 


0.26±0.06 


0.07 ±0.02 


18±6 


19 


WFI J08430294-5302259 


0.63±0.06 


4000 


0.252±0.116 


1.05±0.24 


0.90±0.10 


18±6 


20 


WFIJ08431228-5251411 


3.19±0.19 


2900 


0.002±0.001 


0.18±0.04 


0.07 ±0.02 


64±20 


21 


WFIJ08431860-5306241 


0.18±0.06 


4200 


0.253±0.117 


0.95±0.22 


0.85±0.10 


32±10 


22 


WFI J08433366-5259549 


3.51±0.16 


3000 


0.003±0.001 


0.20±0.05 


0.09±0.02 


64±20 


23 


WFI J08433843-5246379 


2.53±0.06 


3900 


0.076±0.035 


0.61±0.14 


0.60±0.08 


90±40 


24 


WFI J08433939-5307563 


0.02±0.07 


3000 


0.016±0.007 


0.47±0.11 


0.08±0.02 


6±3 


25 


WFI J08433945-5248576 


0.12±0.08 


3300 


0.006±0.003 


0.23±0.05 


0.15±0.04 


80±40 


26 


WFI J08434640-5257036 


0.12±0.10 


4200 


0.180±0.083 


0.80±0.19 


0.75±0.08 


50±20 


27 


WFI J08434751-5250136 


2.00±0.06 


4000 


0.110±0.051 


0.69±0.16 


0.70±0.08 


45±20 


28 


WFI J08434909-5307269 


0.01±0.08 


3200 


0.007±0.003 


0.26±0.06 


0.11 ±0.03 


32±10 


29 


WFI J08440 199-53 10356 


0.41±0.07 


4000 


0.083±0.038 


0.60±0.14 


0.62±0.08 


80±40 


30 


WFI 108404683-5302487* 


0.16±0.17 


4200 


2.128±0.001 


2.77±0.05 


1.40±0.20 


5±2 


31 


WFI J084100I9-5249088* 


0.77* 


3800* 


0.185* 


0.97* 


0.75* 


20* 


32 


WFI 108412176-5247289* 


0.22±0.15 


4300 


1.087±0.001 


1.80±0.01 


1.40±0.20 


6±3 


33 


WFI 108412550-5300154* 


0.13* 


4300* 


0.700* 


2.48* 


1.20* 


10* 


34 


WFI 108413601-5309274* 


0.17±0.07 


3000 


0.012±0.006 


0.39±0.09 


0.07±0.02 


7±3 


35 


WFI 108413888-5259492* 


0.44* 


4100* 


0.101* 


1.13* 


0.65* 


80* 


36 


WFI 108414904-5255192* 


0.05* 


3700* 


0.127* 


0.95* 


0.65* 


16* 


37 


WFI 108415780-5252140* 


0.24±0.2 


4000 


0.203±0.069 


0.74±0.16 


0.85±0.10 


25±8 



■ The uncertainty on Teff is assume d to be 250 K on ave rage (see Sect.l3.lt . 

' From the evolutionary mo dels bv iBaraffe et al. I j 19981) and lChabrier et al I ilOOd) . 

* From the lRobitaille et al.l feood) models for young stellar objects (see Sect. l6.1.Tt . 

* Object showing IR excess emission at 24/im. 

Table 6. Possible multiple systems in IC 2391 (see Sect. [331). 







Primary star 




Separation 






Neighbor star 








ID 


Designation 


Ref.* 


(arcsec) 


R.A. 


Dec. 


V 


Rc 




Previously known 




CI* IC 2391 CTIO 126 


a 


3.1 


08:42:49.385 


-52:52:16.96 


18.25±0.07 


17.13±0.03 


16.46±0.04 


cluster members 




CI* IC 2391 CTIO 25 


a 


2.8 


08:42:46.025 


-52:46:52.64 


20.02±0.07 


18.82±0.03 


17.88±0.06 




6 


WFI 108410814-5255027 


b 


3.9 


08:41:07.726 


-52:55:03.61 


16.00±0.06 


14.98±0.03 


I4.32±0.04 


New member 






b 


3.5 


08:41:08.489 


-52:55:04.19 


15.72±0.05 


15.02±0.04 


I4.70±0.05 


candidates 


26 


WFI 108434640-5257036 


b 


4.5 


08:43:46.814 


-52:57:06.12 


16.39±0.05 


15.68±0.04 


15.26±0.05 




27 


WFI 108434751-5250136 


b 


3.3 


08:43:47.844 


-52:50:12.48 


20.76±0.07 


19.45±0.03 


18.52±0.04 




29 


WFI 108440199-5310356 


b 


2.4 


08:44:02.071 


-53:10:37.92 


15.48±0.06 


14.47±0.03 


13.77±0.04 



t a) lBarrado v Navascues et a Di l2001bh : b) This work. 
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gion of the HR diagram approximately defined by the 10 Myr 
and 100 Myr isochrones (Fig. |5]l- In particular, all our candi- 
dates are confined in the region defined by Tgff =2900-4800 K 
and L*=0. 002-0.4 Lq (i.e. absolute visual magnitude Mv-12- 
6 mag); thus, in order to estimate the contamination level of our 
candidate sample, we have to obtain an estimate of the num- 
bers of stars unrelated to the cluster that may be expected to 
appear in this region of the HR diagram. The number of in- 
terloping stars can be probed by using analytic models of the 
Galactic stellar distribution, i.e. simulations of the expected 
properties of stars seen towards a given direction of the Galaxy 
over a given solid angle. We performed this exercise by us- 
ing the Galaxy model by Robin et al. (2003) and their on-line 
toofl In the temperature range of our candidates, foreground 
stars are expected to be main sequence cool dwarfs, whereas 
red giants are expected to dominate the background popula- 
tion. Assuming a cluster distance of 150-200 pc, we expect 
some 14 foreground dwarfs in the ~30x30 square arcmin area 
observed in IC 2391 with apparent Rc magnitude between 12 
and 21 and spectral types MK, i.e. the ranges estimated for our 
candidates (see Table |4] and |5]). No background giants are ex- 
pected to be found in the locus occupied by the cluster mem- 
bers because all such stars appear much brighter than IC 2391 
members in the same effective temperature range. We thus con- 
clude that only foreground cool main sequence stars can con- 
tribute noticeably to the contamination of our candidate sam- 
ple, with the expected contamination level being 48%. Note 
that this value is very close to the pollution rate of 50% esti- 
mated by lBarrado y Navascues et al. I (l2001bl) for the same area 
in IC 2391. 

Another possible source of contamination is represented by 
galaxies that may have colors similar to those of PMS objects. 
However, such bright galaxies are a priori excluded from our 
selection, because our PSF extraction method allows an effi- 
cient identification/removal of extended objects (Sect. |2.1.4| i. 
The only source of contamination may be the point-like extra- 
galactic obje cts, i.e. QSO's. According to recent QS O's num- 
ber counts ( F ontanot et al J|2007 : Richards et al.ll2005i) . the pre- 
dicted number of QSO's per square degree become significant 
(i.e. greater than ~10/deg^) beyond V ~22 mag, while the ma- 
jority of our candidates are brighter Thus, extra-galactic ob- 
jects might affect our candidate sample by no more than 3%. 

In conclusion, the expected contamination of our candidate 
sample due to field objects is ~50%. 



5. On the sub-stellar population in IC 2391 

The first and only estimate of the IC 2391 mass spec- 
tr um in the very low-mass dom ain has been derived by 



Barrado v Navascues & Stauffej 120031; iBouvier et alj Il998l) . 



Barrado v Navascues et alj (120041) . Independently of the as- 
sumed cluster age between 25 and 50 Myr, they found an index 
of the power-law mass function in the mass range 0.072-0.5 M© 
equal to 0.96±0. 12, consistent with the value reported for other 
open clusters (e.g. iBarrado v Navascues et alj 1200 lal 12002 ; 



^ For details on th e adopted Galaxy model we defer the reader 
to iRobin et"aL I j2003l) . The on-line simulation tool is available at: 
|http://bison.obs-besancon.fr/modele/| 



Below 0.07 Mq, the cluster mass spectrum presents a sudden 
drop below which the authors partially explain by the lack of 
completeness of their survey beyond / a!l8.5 (i.e. 0.05 M© for 
cluster members). 

Although our survey is complete at 100% and 80% lev- 
els down to 0.05 M© and 0.03 M©, respectively, and covers 
60% of the sky area spanned by the cluster, we found only an 
handful of candidate substellar members. In numbers, consid- 
ering the IC 2391 populations updated to January 2008_(i.e. 



IPatten & Simonf 1996 


: Patten & Pavlovski' 1999f |Platais et al. 


2007; Rolleston et al. 


1997; Simon & PattenI 1998), there are 



14 members with Teff <3100K, i.e. below the approximate 
substellar limit for 30-50 Myr old objects. Thus, the fraction 
of sub stellar objects, = yog^JlT^f ' ™ *is clus- 
ter is -8%; by adding our 12 BD candidates (i.e. candi- 
dates with Teff <3100K; see Table |5]), this fraction would 
increase to ~15%. However, all our BD candidates have R- 
magnitudes in the range 15 < R <21, where we expect a 
high level of contamination (70%, see Sect. |4|. Thus, the 
true value of Rss is expected to be very close to 10%. 
Moreover, our physical parametrization indicates that all our 
BD candidates have Tetf >2900 K and masses very close to 
the hydrogen-burning limit. Thus, the results of our deeper 
optical survey confirm the lack of members with spectral 
type later than M7-M8 in IC 2 391 already pointed out by 
iBarrado v Navascues et al.l (l2004t) . A similar behavior has been 
observed in the nearly coeval ope n clusters a Per (~50 Myr) 



and NGC 2547 (-25 Myr) by IBarrado v Navascues et al 



[2002) and lJeffries et al ](l2004. respectiv e ly. Thu s, as already 



concluded by IBarrado v Navascues et al.l (|2004j), one of the 



most solid hypothesis explaining the BD deficit in IC 239 1 and 
other open clusters r emains that by iDobbie et al.l ( |2002|) and 



Jameson et al J (120031) . i.e. the onset of larger-size dust grain for- 
mation in the upper atmosphere of objects with spectral types 
M7-M8 or later, which would be responsible of the local drop 
in the shape of the luminosity-mass relation. 

Of course, a number of alternative hypothesis might be also 
invoked. For instance, simulations of the dynamics of BD pop- 
ulations in open clusters, which consider the effects of a large 
range of primordial binary populations, have also shown that 
the majority of BDs might be contained within primordial bi- 
nary syste ms which then hide s a large proportion of them from 
detection dAdams et al ]|2002h . Alternatively, the deficit of BDs 
in IC 2391 might be due to mass segregation/evaporation ef- 
fects. Dynamical mass segregation acts on a time-scale of the 
order of the relaxation time of a cluster dBonnell & Davies 
1998). which for clusters of the size of IC 2391 is of the 



order of 1 Myr ([Birinev & Tremainelll987l p. 390). Indeed, 
Sagar&Bhat7(ll989h found indication that the higher mass 



stars in IC 2391 have lower velocity dispersions than the lower 
mass stars. 

6. On the debris disk population in IC 2391 

The first s urvey for debris dis ks in IC 2391 was recently per- 
formed bv lSiegler et alj (l2007l) using MIPS on board of Spitzer. 
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The central square degree of IC 2391 was observed at 24/zm, 
with the aim of studying the incidence of debris disks in the 
cluster The authors report IR excess indicative of debris disks 
around 8 previously known members of the cluster, namely 1 
A star, 6 FGK stars and 1 M dwarf, and estimate a debris disk 
fraction of lO^j^ for BA stars and 3l'^l^ for FGK stars. 

In this section we further investigate the debris disk popu- 
lation in IC 2391, using our optic al data in combin ation with 
2MASS data, the MIPS catalog bv'Si egler et"al] 12007) and ad- 
ditional IR data from IRAC. This dataset is suitable to probe 
disk radii from ~0. 1 AU up to 20 AU depending on the stel- 
lar mass (Merm e t al. in preparation). While the work by 
Siegler et alj ( l2007h focuses only on the previously known clus- 
ter members, we extend the search for objects with IR excess 
emission to all the sources detected at 24yum having a coun- 
terpart in our optical WFI images. This allow us to identify 
new candidate KM-type members of the cluster bearing debris 
disks. We remind the reader that our search is spatially limited 
to the overlapping area observed with WFI and Spitzer (Fig.[Tli. 
Moreover, the selected debris disk candidates must be spectro- 
scopically investigated to definitively assess their membership 
to the cluster Indeed, the experienc e with the Spitzer Core to 
Disk Legacy Survey (hereafter c2d; Evans et al.ll2003[) teaches 
us that samples of objects showing mid-IR excess in star form- 
ing regions could be highly contaminated by background post- 
main sequence stars, mostly in the direction of the galactic 
plane. This contamination might account for up to 30% of those 
objects and could be even higher in the case of IC 2391, which 
has a galactic latitude (b«;-7") lower than the typical c2d fields. 
For all these reasons, our search is intended to give a deeper 
insight into the candidate debris disk population in IC2391 . 



6.1. Selection of debris disk candidates 

We identified possible debris disk objects in IC 2391 by mea- 
suring 24fj.m flux densities in excess with respect to the ex- 
pected photospheric em ission; this cr i terion has been adopted 
in the pioneer study by ISiegler et al. I (l2007h and allows us to 
pick-up new debris disk candidates in a homogeneous way. 
We use the {V - Ks) vs. (Ks - [24/ito]) color-color (CC) di- 
agram to identify these excess sourc es across a broad r ange 
of spectral types. As demonstrated by Siegler et al. ( 20071 ). the 
(V - Ks) color is a good tracer for spectral types down to late 
M-dwarfs, while the (Kg - [24fim]) color is a good diagnostic 
for mid-IR excess because the Ks -band flux is essentially pho- 
tospheric for stars older than ~10 Myr, whic h have akeady dis - 
sipated the innermost region of their disks (IHaish et al.ll200ll) . 
We dereddened the magnitudes of all the sources detected at 
24yum using the Ay values derived in Sect. 13. II and the extinc- 
tion law by ICardelh et all ( ll989t Fi We select objects with IR 
excess at 24yum using the mean dereddened photospheric lo- 
cus in_the_(V_- Ks) v s. (Ks - [24jum]) diagram determined 
by ISiegler et aU (|2007|) : for spectral types later than M4, we 
use the I R photospheric colors for M-dwarfs determined by 
Gautier e t al. (2007). 




2 3 
K3-[24/im] 



Fig. 7. (V-ZTs) VS. (/r5-[24;um]) color-color diagram for the 10 
sources in IC 2391 showing IR excess at 24/im. The continuous 
and dashed lines ar e the me an dereddened photospheric loci 
determined by Sieg ler et al.1 12007) and Gautie r et al. (20o3), 
respectively. Sources redder than this threshold possess (Ks - 
[24/im]) flux ratios in excess of expected photospheric colors. 
Symbols are as follows: the stars mark the confirmed cluster 
members, while squares and crosses mark the member can- 
didates and the four sources of dubious nature discussed in 
Sect. 16.1.11 respectively. 



As shown in Fig. |7] we found 10 sources redder than this 
threshold, i.e. objects possessing (Ks - [24yum]) flux ratios in 
excess of expected photospheric colors by at least 15%, i.e. the 
same criterion as in Siegler et al. (2007). Four out of these 10 
sources (ID 30, ID 32, ID 34 and ID 37), are previously known 
members of the clusteiQ (stars in Fig. |7]i, two (ID 1 and ID 7) 
are new candidate members selected in this work (squares in 
Fig. |7|i, while the remaining four (ID 31, ID 33, ID 35 and 
ID 36) are new sources with IR excess emission neither re- 
ported by previous studies nor selected as candidate cluster 
members by our optical criteria (crosses in Fig. |7]i; the latter 
objects are further discussed in Sect. I6.1.TI 

The optical and near-IR photometry and the physical pa- 
rameters of the 10 sources showing IR excess at MIPS/24/im 
are reported in Tableland Table|5] respectively, while Table|2] 
reports their IRAC and MIPS/24//m magnitudes. 

6.1 .1 . Comments on individual objects 



We adopt the lCardelli et all ( Il989l) extinction law for consistency 
with our parametrization procedure (Sect. lJTt . 



Siegler et al.l ( 120071) reported IR excesses at lA-jim indicative of 
debris disks for 8 objects in IC 2391. Only two of these objects 
fall in the area observed with WFI, namely ID 29 and ID 24 
from their Table 2; we correctly recovered ID 29 as a possible 
debris disk object (i.e. our ID 37), while ID 24 is saturated 
in our WFI images and was then a priori excluded from our 
analysis. 

Four of our debris disk candidates, namely ID 31, ID 33, 
ID 35 and ID 36, though showing a clear IR excess emission 
in their SEDs (Fig. [8]l, were not selected as candidate cluster 



See Table|4]for previous bibliography on these objects. 
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Table 7. IRAC and MIPS/24yum magnitudes and SED slope ia[K&24fim]) for the 10 objects showing IR excess emission at 24//m 
identified in this work. 



ID 


Designation 


Status'' 


lRAC/3.6yum 


lRAC/4.5//m 


lRAC/5.8yum 


lRAC/8.0/im 


MlPS/24yum 




1 


WFl J08403964-5250429 


HR-C 


14.26+0.01 


14.02+0.01 


14.13+0.08 


14.11+0.13 


11.07+0.25 


-1.63 


30 


WFl J08404683-5 302487 


M 


15.76±0.04 


15.61+0.06 


15.92+0.34 


15.90+0.98 


10.71+0.18 


-0.89 


31 


WFl J08410019-5249088 


IR-C 


12.46+0.01 


12.54+0.01 


12.47+0.02 


12.39+0.04 


10.31+0.14 


-2.07 


7 


WFl J08410886-5248387 


HR-C,X 


12.32±0.01 


12.21+0.01 


12.22+0.02 


12.10+0.03 


11.42+0.47 


-2.45 


32 


WFl J08412 176-5247289 


M 


15.58±0.03 


15.59+0.05 


15.26+0.23 


16.88+1.60 


11.41+0.32 


-1.44 


33 


WFl 108412550-5300154+ 


IR-C 


15.88+0.05 


15.93+0.06 






10.54+0.17 


-0.96 


34 


WFl J08413601-5309274 


M 


12.35+0.01 


12.24+0.01 


12.14+0.02 


12.17+0.03 


11.44+0.30 


-2.46 


35 


WFl J08413888-5259492 


IR-C 


13.70+0.01 


13.58+0.02 


13.28+0.05 


12.95+0.06 


9.86+0.09 


-1.34 


36 


WFl J08414904-5255192 


IR-C 


14.07+0.02 


13.83+0.02 


13.44+0.06 


13.03+0.06 


10.20+0.15 


-1.26 


37 


WFl J08415780-5252140 


M 


10.18+0.01 


10.22+0.01 


10.17+0.01 


10.10+0.01 


9.73+0.15 


-2.71 



^ M = confirmed cluster member; HR-C = member candidate according to the location on the HR diagram; X = X-ray sources iMarino et alj 
I2OO5I) ; IR-C: young stellar object candidate according to the IR properties (Sect.[6jT]l. 
+ Flux at 5.8 and 8 fim not available. 



members neither by previous surveys nor by our optical cri- 
teria. This is because these four sources are sub-luminous in 
the HR diagram when adopting the distance of 150 pc and the 
extinction determined as explained in Sect. 13.11 Possible rea- 
sons for this inconsistency may be an underestimate of the ex- 
tinction and/or the distance for these four sources. From a vi- 
sual inspection of our optical WFl images, these four sources 
appear point-like and are concentrated in the central area of 



the cluster (Fi 



RobitaiUe et 



tng. 
aU 



Tj); ni oreover, by fitting their SEDs with the 
( 2006r) grid of models for young stellar ob- 



jects normalized to the observed /-band flux, we obtain stellar 
and disk parameters still consistent with young stars belonging 
to the cluster (see Table |5]l. Thus, for the time being, we in- 
clude these sources in our list of possible debris disk objects in 
IC 2391; future follow-up spectroscopy will shed light on their 
nature. 



6.2. On the disk properties 

This section aims for a description of the circumstellar ma- 
terial around the 10 objects with IR excess emission selected 
above, in order to give some hint on the nature a nd evolution 



ary phase of their disks. Recent pubhcations ("e.g. lAlcala et al 



20081; iHarvev et al.ll2007l; iMerfnet al.lboO S') of Spitzer obser- 
vations in star forming regions within the frame of the c2d 
survey have provided several methods to study disk properties 
such as the amount of circumstellar dust, morphology, etc. We 
applied some of these methods to further investigate the nature 
of our 10 debris disk candidates in IC 2391. 

We first construct the SEDs of these objects from optical to 
mid-IR wavelengths, as aheady described in Sect. I3.1l but in- 
cluding IRAC and MIPS data and, wher e possible, additional B 
magnitudes from the NOMAD catalog ( Zacharias et alJE005 ). 
For each object the flux at each wavelength was corrected for 
interstellar extinction using the Ay values reported in Table |5] 



and th e normal (Rv=3.1) reddening law by ICardelli et al. 
19891 ). The SEDs of the 10 debris disk candidates in IC 2391 



are sh own in Fig.[8ltogether with the best-fittmg mo del spec- 
tra bv lHauschildt et all d 19991) and lAllard et"ai[| (I2OOOI) with the 
same temperature as the objects (see Sect. 13.1b . representing 
the stellar contribution. For the majority of the objects the IR 
excess with respect to the expected photospheric emission is 
clearly visible. Moreover, the lack of near-IR excess emission 
suggests the presence of disks whose inner part has been al- 
ready dissipated (see, e.g. JSiciha-Aguilar et alJE006l) . 

Additional clues on the nature of our debris disk candi- 
dates can be obtained by investigating the o bjects IR classes 
according to the de finition by ,Ladal ( 1987*) as extended by 
Greene et al.l (119941) . This classification is traditionally based 
on the value of the a spectral index, which is determined as 
the slope of the SED in the form log AF^ vs. log A at wave- 
lengths longer tha i i 2tim. Following the c2d convention (e.g. 



Alcala et al. 2008; Harvey et al. 2007; Merfn et al. 2008), we 



use a[K&24iAm] as the reference spectral index, which is the slope 
of the linear fit to the fluxes at Ks and MIPS/24yL(m (Table|7]i. As 
shown in Fig.|9] we found that 5 out of the 10 debris disk candi- 
dates in our sample present a[K&24jjm] slopes typical of IR class 
III sources (aiK&24u m} <-1.6) according to the classification by 
Greene et al.l ( 119941) . i.e. they have optically thin disks. The re- 
maining 5 objects fall in the class II regime (- 1 .6< a[K&24nm] ^- 
0.3), i.e. they would have optically thick disks; however, the 
majority of them lay very close to the class III boundary (i.e. 
a[K&24i.im] ^-1)- The possibility that IC 2391 would have any 
thick/transitional disks is quite surprising considering its age 
(30-50 Myr), because these disks, along wit h primordial disks. 



are typically dissipated within ~10 Myr (IDecin et al.l 12003 



Dominik & Decinll2003h . However, Spitzer data have shown a 
very clear mass dependence of the disk dissipation timescale, 
i.e. very low-mass stars dissipate th eir primo rdial disks slowe r 
than more massive stars ([Hillenbrand 2008; Scholz et al ■l2007h . 
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Fig. 8. Spectral energy distributions of the 10 sources in IC 2391 showing IR excess emission identified in this work. The 
dereddened fluxes are represented with filled circles, while the obs erved fluxes are represe nted with open circles. Flux upper 
limits are marked with an arrow. Th e best-fit NextGen spectrum by 'H auschildt et aL ( 199 9*). for objects with Teff >4000 K, or 
StarDusty spectrum by Alla rd et al. I IJOOO), for objects with Ten <4000 K, is overplotted on each SED, representing the stellar 
emission. Each panel is labeled with the object name (and ID number between brackets), estimated effective temperature (Table|5]l 
and membership status (see Table|7]). 



In conclusion, the 10 sources with IR excess emission se- 
lected in Sect. 16.11 are likely to possess evolved disks, whose 
inner part might have already been cleared by the formation of 
large-sized bodies (such as rocks). 

7. Conclusions 

We performed a deep wide-field imaging survey of the young 
open cluster IC 2391 to study its low mass population down 
to 0.03 Mq. We used our optical VRI data, in combination 
with JHKs magnitudes from the 2MASS catalog to estimate 
the physical parameters of the optically detected sources and 
selected cluster member candidates on the basis of their loca- 
tion on the HR diagram compared to theoretical isochrones. We 



estimated the contamination using statistical arguments and, 
when possible, proper motion from the NOMAD catalog. Our 
survey has a completeness level of 80% down to /c=23.0, i.e. 
~l-2 mag fainter than previous optical surveys conducted in 
IC 2391, and revealed 29 new low-mass member candidates, 
among which 12 may be BDs. The expected contaminations is 
at least -50%. 

We confirm the presence of a significant drop of the cluster 
mass spectrum across the stellar/substellar boundary, with the 
fraction of sub-stellar objects being of only 8-15%. The lack 
of BDs in IC 2391 might be due to the dip in the luminosity 
function around M7-M8 already observed in other star-forming 
regions, open clusters and the field, which is caused by the be- 
ginning of larger-size dust formation in the atmospheres of ob- 
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Fig. 9. Distribution of the a[K&24nm] SED slopes for the sources 
in IC 2391 showing IR excess emission identified in this work. 
The dashed vertical lines mark the boundaries of class I, flat, 
class II and class III objects according to the cl assification by 
Ladal ( Il987l) as extended bv lOreene etal] (Il994 



jects in this efifective temperature regime. However, other ex- 
planations such as dynamical mass segregation and/or masking 
in binary systems might be invoked. Our data are not sufficient 
to decide between these scenarios. 

Finally, we have combined our ground-based optical ob- 
servations with near-IR photometry from 2MASS and Spitzer 
GTO observations to investigate the debris disk population of 
IC 2391. We identified 10 possible debris disk objects in the 
cluster on the basis of their lAjjm flux densities in excess with 
respect to the expected photospheric emission. We constructed 
optical/IR SEDs for all these objects and investigate the prop- 
erties of their circumstellar material. According to our analysis, 
these 10 objects are likely to possess evolved disks whose inner 
part has been akeady cleared up. A few of these objects might 
possess thick/transition disks, what would be surprising, given 
flieage ofIC 2391. 
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